Singular response to a dopant of an evaporating crystal surface 
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Moving crystal surfaces can undergo step-bunching instabilities, when subject to an electric cur- 
rent. We show analytically that an infinitesimal quantity of a dopant may invert the stability, 
whatever the sign of the current. Our study is relevant for experimental results [S. S. Kosolobov et 
al, JETP Lett. 81, 117 (2005)] on an evaporating Si(lll) surface, which show a singular response 
to Au doping, whose density distribution is related to inhomogeneous Si diffusion. 
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Introduction. — Crystal surfaces may undergo dynam- 
ical instabilities while growing or evaporating^ 2 .. An 
important class of instabilities is that determined by 
an asymmetric current of adatoms diffusing on terraces. 
This asymmetry may be intrinsic, e.g. determined by 
asymmetric attachment to steps 3 - or by surface recon- 
struction^, or it may be extrinsic, e.g. due to impurities^ 
or to an electric current^. Extrinsic instabilities may have 
the advantage to be tunable and therefore to be suitably 
switched on and off. In this Letter we consider an ex- 
ample of extrinsic instability, displaying a phenomenon 
which is of double interest. An interest for fundamen- 
tal noncquilibrium physics, because we propose a model 
showing a singular response to doping. But also an in- 
terest for applications, because we suggest that recent 
experimental results^ on evaporating Si(lll) surfaces un- 
der an electric current, showing a singular response when 
exposed to a variable quantity of Au, can be understood 
at the light of our model. 

Vicinal surfaces are obtained cutting a crystal along an 
orientation which is close to a high symmetry one, result- 
ing in a morphology similar to a flight of steps. During 
evaporation, atoms detach from steps and diffuse on ter- 
races until possible desorption. The application of an 
electric current induces an elcctromigration force, whose 
direction may depend on temperature T: this is the case 
for the widely studied case of Si(lll), where a sequence 
of stable and unstable regions are found when varying 
T— . Here we are interested to analyze how stability is af- 
fected by doping. We assume (and later we discuss) that 
doping induces inhomogeneity in the diffusion process, 
both in its symmetric part and in the drift term, the lat- 
ter modifying the elcctromigration force. In particular, 
we show that a divergent drift in a region of vanishing 
size may induce a change of stability, even if the sign of 
the drift does not change. This result might not appear 
surprising as it is, because step decoration^, changing at- 
tachment kinetic coefficients may have a similar effect. 
However, step decoration misses two important features: 
first, in our model and in the Si(lll) experimental sys- 
tem, doping induces a change of stability for both signs 
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FIG. 1. (a) Profile of a one dimensional surface and (b) pro- 
file of a terrace, relevant for the inhomogeneous model. All 
symbols are explained in the text. 



of the current, meaning that doping cannot be reduced 
to a fixed stabilizing or destabilizing effect; second, an 
increasing quantity of dopant may induce further stabil- 
ity/instability transitions. In this Letter we focus on the 
very first transition, occurring at a vanishing critical den- 
sity of the dopant. However, we also argue how a larger 
quantity of dopant may further change the stability of 
the surface. 

The model. — The simplest model of adatom diffusion 
in the presence of evaporation and drift (see Fig. QJi) re- 
quires to solve the following stationary diffusion equation 
for the adatom density p(x), 



Dp"(x) - 2a^ 1 5p'{x) - i 2 p{x) = 0, 



(1) 



where ao is the lattice constant; D is the diffusion con- 
stant; 8/ao = q a{{ ED / (ksT) is the drift, with q ett being 
the effective charge of adatoms and q g{f E being the elec- 
tromigration force; 7 2 = 1/t is the evaporation rate. As 
illustrated in Fig. [T|d, we adopt a simple model for inho- 
mogeneity: each terrace is separated in two regions with 
different values for D and S. Boundary conditions are 
expressed in terms of the surface current 



J(x) = -Dp'(x) + 2ao 1 Sp{x), 



(2) 
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which at steps is proportional to the supersaturation, 
J stcp = ±K(p atop - p cq ). 



(3) 



Here the plus (minus) sign applies to the ascending (de- 
scending) step. We expressly chose the same kinetic coef- 
ficient K for ascending (K + ) and descending (K_) step, 
to avoid stabilizing or destabilizing effects due to their 
asymmetry, K + 7^ K _— 

The general solution of Eq. ([I) for the piece of terrace 
x € [a£, £] is given by 



p(x) = Ae 
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and analogously for x € [0,a£], with coefficients A' , B' 
in The unknowns A, £?, A' , B' are determined from 
two conditions (j3J) at steps x = 0, £ and from continuity 
of current ([2]) and density p(x) at a; = a£. Before going 
further, let us remind what is the relevant function de- 
termining the stability of the evaporating surface. Each 
step moves with a velocity proportional to the sum of the 
upper and lower step current, J stop , as given by Eq. (|3|). 
Perturbing a perfect train of steps, we findi£ the stability 
being governed by the function (j){£) = p(£) — j»(0): the 
system is stable (unstable), if de(f>(£) is positive (nega- 
tive). 

We now assume a diverging drift in the region close 
to the descending step (5' — > 00), whose size vanishes in 
such a way that the product 



S'a£ 
a D 



(G) 



is constant. This model, as discussed at length in the 
experimental Section, is motivated by a possible strong 
impact on drift of a dopant higly localized near steps. 
As shown in detail in the Supplemental Material^, we 
don't really need a "singular" model with a diverging 5' 
in order to get reversal stability, but the experimental 
system does show a singular behavior. 

We have also considered the case of a diverging D, 
but it results to have no effect when restricted to a re- 
gion of vanishing size. For this reason, the parameter D' 
does not appear in the following equations. For ease of 
notation, we define the scaled quantities 5 = <5/(ao£>), 
7 = "f/V~D, and K = K/D, all of dimension [X -1 ]. 

Solving the two coupled diffusion problems in the 
above limit, we have found that the resulting effect is 
equivalent to introduce a new boundary condition at the 
descending step (x = 0), 
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The new condition arises because the density profile in 
proximity of the descending step becomes infinitely steep 



in the limit a£ — ¥ 0, changing from some value p 5top at x = 
to a value p Btcp e 2c at x = a£ — > + . Therefore, p stop e 2c 
becomes our new postulated value p(x = 0), after taking 
the limit a£ — > 0. Because of this redefinition of J st c P , 
the function (/>(£), which governs the stability of surface 
dynamics, has to be redefined as 4>(£) = p(£) — p(0)e~ 2c . 

By solving Eq. (jTJ) with usual boundary condition ([3]) 
at x = £ and new boundary condition ([7]) at x = 0, and 
assuming 



we obtain 



) 



«<7^< 1, 



2K 2 e- 2c p ccl (4:S - K(l - e~ 2c )) 



(8) 



(K + Ke- 2c + 2SK£ + K 2 e~ 2c £) 2 
Therefore, the instability condition de<fi(£) < reads 



(9) 



45 < K(l 
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(10) 



which reduces to the usual condition 6 < 0, when bound- 
ary layer is absent (c = 0). It is worthnoting that Eq. (TTU)) 
does not depend on the desorption rate 7. In fact, it is 
possible to study a model where evaporation is neglected: 
in this case, simpler calculations allow to keep 5' large, 
but finite. This calculation, which gives the same result 
and therefore proves the robustness of our model, can be 
found in the Supplemental Material^. In the experimen- 
tal system^ we consider here, desorption is not negligible, 
see Eq. ([5]) , but we have shown it does not affect the sta- 
bility reversal process. 

If Eq. ([TO"]) is satisfied, the instability time, T inB , is set 
by the relation^ 



a K\diqb{£)\. 



(11) 



Analogously, one may consider a situation with di- 
verging drift in the region close to the ascending step, 
8'£(1 — a)/(aoD) = ci as a — > 1. This inhomogeneity 
results in a new effective boundary condition at the as- 
cending step, 



J.tcp = K(p stcp e 



2ci 



(12) 



and in the redefinition <f>(£) = p(£)e 2ci —p(0), so that the 
instability condition dt4>{£) < becomes 



45 < K(e 2ci - 1) 



(13) 



Expressions (jTDJ) and (|T5)) arc a central result of our 
paper. In the following we discuss their application to a 
surprising experimental effect, observed on the evaporat- 
ing surface of Si(lll) in the presence of a tiny coverage 
on Au atoms. 

The experiment^. — Step bunching instability of a 
Si(lll) vicinal surface heated by a DC electric current to 
high temperatures is a well established phenomeno n 13 ' 14 . 
In the so-called first temperature regime (830°C< T < 
950°C), the regular vicinal surface is stable (unstable) 
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for an uphill (downhill) current. The effect is usually 
understood^ by assuming that neutral silicon adatoms 
acquire an effective positive charge which gives rise to an 
electromigration force S. In the first temperature regime, 
approximations ((5J are valicU^, the density profile is lin- 
ear, p(x) ~ 5x, and we have stability for S > 0. 

A few years ago, it was observed by Kosolobov, Laty- 
shev and collaborators (KL) that a submonolaycr deposi- 
tion of gold on a vicinal Si(lll) surface drastically affects 
its stability, which changes four times as a function of the 
increasing Au cover agei. In particular, the addition of a 
very small quantity of Au atoms to a clean Si surface 
(0.0016 ML in RefJ£) resulted in a reversal of stability 
(uphill DC setup becomes unstable and viceversa). Such 
critical density, not enough even to decorate vicinal steps 
with a single line of Au atoms, hints at a phase transition 
at zero critical density, p A T u — > as I — > oo. Moreover, 
the instability develops much faster than is usual at this 
temperature for a clean Si surface^. 

We propose to relate this phenomenon to our model of 
Fig. 1(b), with a diverging drift. To do so, we note that 
Au adatoms have a much larger effective charge than 
Si adatom o 18 i 19 , so that their dynamics in the electric 
field is faster. As a consequence, Si adatoms move in an 
established environment of Au adatoms, whose steady 
profile may be strongly inhomogeneous near the border. 
Relating Au inhomogeneity to Si diffusion inhomogeneity 
causes a stability reversal, as demonstrated below. 

The key observation is that the negative effective 
charge of Au atoms and their strong affinity^ to the step 
region, make them concentrate near the descending step 
boundary x = 0. The effective equilibrium profile of Au 
adatoms, n(x), is governed by an equation of type ([l), 
i.e. 
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FIG. 2. Schematic profiles of the density of dopant adatoms 
(thick lines) and drifting forces S(x) (thin lines) at small 
dopant concentration, p Axi l -C 1. For uphill (downhill) cur- 
rent S > 0( S < 0), the dopants are concentrated near the 
left (right) boundary. The S < curves are shifted for better 



near the upper step boundary, 
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As anticipated, we relate the gradient in the density of 
Au atoms to an additional drift (effective electromigra- 
tion force) of Si adatoms as — 



S' = S- 



„ dn{x) 
dx 



(18) 



D Aa n"(x) + 2a,7 1 | ( 5 Au |n / (a;) = 0. 



(14) 



Due to strong affinity assumption, the boundary con- 



dition is n(0) 



where nj u is the equilibrium density 



of Au adatoms close to the descending step. Within the 
experimental instability time scale, Au coverage remains 
approximately constant. Therefore, we can solve (|14p 
assuming a negligible Au evaporation and get n(x) = 
n+e" KI , where k = 2a ~ 1 \5 All \/ D Au . 

If we now impose that the space integral of n(x) equals 
the total number of deposited Au atoms, 



Pau^ = / n{x)dx = — —, 
Jo « 



(15) 



we get the relation n = n^ u /p Au £, which diverges for 
Pau — > 0. Finally, the equilibrium density profile writes 

n(z)=n+ expf-^jf). (16) 

Analogously, for downhill current and small density 
of Au adatoms we obtain a distribution function peaked 



This rcnormalizes the drift in the region where the gra- 
dient is sufficiently large (see Fig. 2), getting 



S(x) = 5' = c*n Avi K, for x G [0, K 



(19) 



If we compare Eq. (fTi?|) with Eq. © and identify the size 
(a£) of the vanishing region of divergent drift with K , 
we find the relation c = c*(n+ u /ao-D). In conclusion, we 
have related the problem of Si diffusion in the presence 
of a tiny quantity of Au doping to the model of inhomo- 
geneous diffusion, discussed in the previous part of this 
Letter. 

We stress that our results do not change qualitatively, 
if the parameter c has a slow dependence on other quan- 
tities such as the local gradient or average density of Au 
adatoms: the decisive role is played by the sign of c. Note 
also that the limit of a diverging drift and gradient is a 
simplifying assumption chosen for clarity of presentation. 
Simulation of a discrete model with strictly finite drifts 
and gradients yields qualitatively the same results as the 
continuous model (to be discussed elsewhere). 

Discussion. — The parameters given in Refs i 22 i 23 lead 
to dimensionless estimates 8£ w 10~ 6 , ; y£ w 10~ 4 and 
K£ sa 10 -2 , where £ ps 100 nm is a typical size of a 
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terrace. This means that assumptions (JSJ are satisfied 
and Eqs. (| MT5)) do apply to the KL experiments. 

For a tiny quantity of gold, Eqs. (|TU|) and (|TT|) should 
be evaluated. Since K/8 ps 10 , Eq. (fT0|) is surely satis- 
fied for finite c (more precisely, for c > 2 x 10 -4 ). The 
microscopic origin of the strong Au-induced increasing of 
the effective drift S(x) for silicon adatoms may be a hard 
core exclusion interaction between Au and Si adatoms: 
the presence of Au adatoms near the step excludes the 
presence there of Si adatoms and thus generates an ef- 
fective enhanced drift of Si adatoms close to the step2i, 
proportional to the local density gradient, as in (|18[). 

Another microscopic origin might be the recharging ef- 
fect 1 ^, according to which a neutral Au adatom subtracts 
negative charge from the substrate, while the substrate, 
to compensate, transfers equal charge of the opposite sign 
to the Si adatoms. Since qgf w 0.004e±£ and qf u « -e±2, 
we expect that recharging effect may lead to a strong 
renormalization of 8. Note that the drift enhancement 
due to recharging is proportional to the local density of 
the Au adatoms itself, and not to its gradient as in case of 
the hard-core exclusion effect alone. With some amend- 
ments to our line of argument (fl~5j) - (flT))) . the instability 
reversal can be obtained. We note that both recharging 
and hard core exclusion effects contribute with the same 
sign to the renormalization (TTS)) , leading to the discussed 
above kinetic instability reversal. We also stress that the 
effective charge of Si adatoms remains positive across all 
the transitions. 

A piece of evidence that c is not small and also that our 
theory does apply to KL experiments is the expectation 
that Au doping leads to a much stronger instability than 
simple reversing of the sign of the current. In fact, since 
K(l - e~ 2c ) > 5, according to Eqs. (|9I11|) we have that 
T do P i„g ^ T rovcrsai^ ag ^ j g actually seen in experiments 1 ^. 

In the Introduction we have stressed that the same 
doping induces a change of stability for both signs of 
the current. Therefore, let's now consider a downhill 
current, which means 8 < and that the clean Si sur- 
face is bunching-unstable, see Eq. (flU)) with c\ = 0. 
Au adatoms are now driven towards the ascending step 
and the boundary layer (inhomogeneity) forms at the as- 
cending step, with a profile given by Eq. (|17p (see also 
Fig. [5]). Assuming an additional drift of Si adatoms, 
8' = —c* = c*n~ u K, in the region of a strong Au gra- 
dient, (l — x) < 1/k, we recover our model with diverging 
drift 8' at the ascending step. The same considerations 
we did for positive 8 and Eq. (fT0| are now applicable to 
negative 8 and Eq. (|13j). providing perfectly symmetrical 
conclusions for a negative electromigration force: a tiny 
amount of Au doping allows to stabilize the surface and 
the stability is much stronger than the stability gained 



by simply reversing the field. 

Let's now pass to reason on further features of KL ex- 
periments^, which go beyond our simple model. These 
experiments are done for a variable quantity of Au, from 
a clean Si surface (p Au = 0) to an almost full covering 
{Pau = 1)^ and authors find a total of four stability 

1 ^ u 

FIG. 3. Stability diagram of a Si(lll) surface, with an up- 
hill current (S > 0), when doped with a variable quantity of 
Au. The clean surface (p Au = 0) and the fully covered Si 
surface (p Au = 1) are stable. A tiny Au deposition on the 
clean surface or a tiny removal of it from the fully covered 
surface induces a stability reversal and the surface becomes 
unstable (thick horizontal lines). In the intermediate region, 
the surface is expected to recover stability (thin horizontal 
line), because the Au profile is more flat. If the sign of the 
current is reversed, stability regions are also reversed (thick 
line thin line). 

transitions, while increasing p Au (see Fig. 3). It is rea- 
sonable to assume that enhanced drift is suppressed with 
increasing Au coverage, which becomes more homoge- 
neous. Since a model with a constant 8(x) is equivalent 
to a clean Si surface, we expect that increasing p Au leads 
to a new reversal of stability, in agreement with experi- 
ments. However, why a further increase of /9 Au produces 
a further change of stability? This change of stability 
is better understood if large p Au coverages are described 
in terms of "Au holes" , i.e. empty places which can po- 
tentially host Au adatoms: p A ° u lcs = 1 — p Au . A gradient 
of Au adatoms corresponds to a gradient of Au holes, 
which appear to have a strong positive effective charge. 
This change of sign compensates the minus relating the 
two spatial derivatives, /3 A ° los = —p' Aa - Thus, removing 
a tiny quantity of Au from a fully covered Si surface re- 
verses its stability. In Fig. 3 we graphically summarize 
the stability diagram of Au-doped Si(lll). 

We conclude with one prediction of our model: if 
dopant adatoms have a positive effective charge, such as 
Cu or Agi£, the stability reversal effect is not expected. 
Indeed, positively charged foreign adatoms under uphill 
current (8 > 0) will be driven to a ascending step and 
eventually form a profile (jTTJ) , for which (fill)) will apply. 
Proceeding along the lines (|18til9[) , we conclude that (|T3|) 
cannot be satisfied, since 8 > and c\ < 0, so that an 
instability reversal will not happen. Indeed, for a Cu- 
doped surface, the reversal of bunching stability was not 
observed 1 ^. 
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